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CN ' We thank Park et al. for their stimulating comments |1|]. 

Oj, Notice the following points (in the whole Reply we will adopt the same notation as in |l|,0l): 

D ' i) It is important to observe that Hi is not only quartic, but involves also a quadratic term {L'^ (l — cosq)aqa^q, 

^^ , in their notation) which is essential for the probability conservation. It seems that Park et al. have forgotten this 

^\ ' contribution in their computations. 

^~^ ii) The expression (2) of their comment |l| is incorrect and should read 

I— I Z[^,tQ + dt] = exp I ^ ^Q^g^_q[/(q,io) — ciil'^l'Z) + ^*(9) — 2esinq))/((7,io) + 2e(iisinq] j , where we use the same 

'~^ notation as Park et al. 

fl^ . iii)Even if we have some doubts on the correctness of the explicit form of their generating functions Z[^, di], we 

Ch ' agree that in the general case the latter involves a non-gaussian form. Our original approach follows, however, a 

I . different path. 

jrt ' Let us illustrate, in the simplest case, what we meant in writting Ig]: "Because of the conservation of probability, 
an analog of the Wick theorem applies and all multipoint correlation functions can be computed" . To do so we 
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recall that H = Hq + j^Hi, where Hq = J2n>Q '^(9)flJog + u!*{q)aLga^q + 2smq 



"q^-q I ^ \HJ"'-q^-q < '^ "'^^ y \^ ^qO-^q + €.a_qaq 



eL and 



Hi = -LJ2qi^ - cosq)a\aq + J^q^.q^.qa cos{qi - q2)a\^aq^a\^aq^+q^_q^ , where Hq is the free-fermionic part of H. 
Now, because of the conservation of probability, we have {x\H = 0. Morever, Hq being itself stochastic, we also have 
^ ; {x\Ho = {x\Hi = 0. 
^ ■ This implies then, 

o ; 

,^ '. L'^il- cos q){x\alaq^ ^ cos{qi - q2){x\al^aq^al^aq,+q^-q^ (1) 

9 qi,q2,q3 

►^ _ Because of the translational invariance {q ^ 0), we have {x\a}jaq'e^^*'\p) — {a^qaqi){t)5q^qi — g{q,t). We also define 

\0' F4t)^j:^^giq,t)cosnq. 

C^ ! It follows from (0), that 

(N. ^ 

^; l2 XI ^os{qi -q2){alfy^al^aq,+q.,-q^){t) = J- J2{1~ cos q){alaq){t)=Foit) -Flit) (2) 

f^ ' 91,12, qa q 

Now, the translational invariance implies that 



(n™n„,+i}(i) = — Y^ cos{qi - q2){al^aq^al^aq,+q^-q^){t) (3) 



91,92,93 

We then obtain the following expression for Ci{t) (with the same notation as in B): 



X3 

O ■ Ci{t) = {n„,n„,+i){t) - {p{t)f = Fo(i) - Fi{t) - {Fo{t)f (4) 



Such an expression, which only comes from the probability conservation, can also be obtained via an "a la WicK^ 
factorization of the correlation functions, as given by Eq.(9) of Ref. |0], and this independently of the form of Z[^, t]. 
Let us close this Reply with the following remarks: 



j^ ' iv) Our approach provides, in the model considered, results which are certainly beyond traditional Hartree-Fock 

theory. In fact the latter predicts decay of the density as i"'^, for the model under consideration (when e = 0), instead 
of the correct behavior, pit) ~ (Airih + h')t)~^^^ — ^{Tre'{h + h')t)^^ |3|, which we can reproduce (it has to been 
emphasized that we obtain not only the leading term of p{t) but also the subdominant one), 
v) V7, we are able to provide the exact and non-trivial steady-states Q. 
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